Accepted for publication in ApJ 

Preprint typeset using I^T^X style emulateapj v. 03/07/07 



SPATIALLY RESOLVED SPECTROSCOPIC STAR FORMATION HISTORIES OF NEARBY DISKS: HINTS OF 

STELLAR MIGRATION* 

Peter Yoachim\ Rok Roskar^, Victor P. Debattista^ 

Accepted for publication in ApJ, March 1, 2013 

ABSTRACT 

We use the Mitchell Spectrograph (formerly VIRUS-P) to observe 12 nearby disk galaxies. We 
successfully measure ages in the outer disk in six systems. In three cases (NGC 2684, NGC 6155, and 
NGC 7437), we find that a downward break in the disk surface brightness profile corresponds with a 
change in the dominant stellar population with the interior being dominated by active star formation 
and the exterior having older stellar populations that are best-fit with star formation histories that 
decline with time. The observed increase in average stellar ages beyond a profile break is similar to 
theoretical models that predict surface brightness breaks are caused by stellar migration, with the 
outer disk being populated from scattered old interior stars. 

In three more cases (IC 1132, NGC 4904, and NGC 6691), we find no significant change in the 
stellar population as one crosses the break radius. In these galaxies, both the inner and outer disks 
are dominated by active star formation and younger stellar populations. 

While radial migration can contribute to the stellar populations beyond the break, it appears more 
than one mechanism is required to explain all of our observed stellar profile breaks. 
Subject headings: galaxies: kinematics and dynamics — galaxies: formation — galaxies: structure 



1. INTRODUCTION 

The stellar components of spiral galaxies are typically 
found embedded in much larger HI disks (Bosma 1981; 
Broeils & Rhee||1997| [Begum et al. 2005) [Walter era | 



2008p. Early work based on photograph ic plates found 
disk galaxies to be sharply truncated ([van der Kruit 
19791 Ivan der Kruit fc Searle[[1981b|a[ ). Recent work by 
Pohlen Trujillo (200lf based on Sloan Digital Sky Sur- 
vey (SDSS) images finds three main types of disk galaxy 
surface brightness profiles. The majority of galaxies in 
their sample (60%) are best fit with an exponential pro- 
file that transitions to a steeper exponential (at what is 
commonly referred to as a break radius). Another 30% of 
galaxies have profiles which transition to shallower scale- 
lengths (up-bending profiles). Only 10% of galaxies are 
well fit by a single exponential profile. 

The first step in understanding disk outskirts therefore 
needs to be an explanation for the down-bending profiles. 
Traditionally, it has been assumed there is a surface den- 
sity thresh old below which gas will fail to undergo star 
formation (Kennicutt 1989 Schaye 2004). There is a 



long hist oryof hunting for star formation at large galac 



tic radii. Ferguson et al. ( 1998 ) discovered extended star 
formation in 3 galaxies based on Ha imaging. Large sam- 
ples of galaxies have been imaged with GALEX, show- 
ing that 20% have UV-bright emission features beyond 



the traditional star formation threshold r adius ( Thilker 
eTaT] [20051 [Zaritsky fc Christlein| [20071 [l^hilker et aL 
2007D. One interpretation of the outer disks in these 
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systems is that, albeit at lower effi ciency, th e majority 
of the stars formed lo cally (Elmegreen fc Hunter|^2006} 
However ^' ' ' 



Bigiel et al. (2010) find that star formation is 



extremely inetticient m the extended HI dominated en- 
velope. Thus the majority of stars beyond the break in 
down-bending profiles probably did not form in situ. 

Ros kar et al.^ (2008b ) present simulations that show 
how the combination of a star formation threshold and 
secular evolution can create disks with broken exponen- 
tial profiles. In this scenario, the region beyond the star 
formation cut-off is populated by stars which migrate 
there from the inner disk. Because migration is a ran- 
dom process, with stars moving both outwards and in- 
wards, it takes stars progressively longer to reach larger 
radii. Thus this process creates an increasing mean stel- 
lar age beyond the break radius. Even if there is some 
low-level star formation in the outer disk, the migrated 
popul ation can still be expected to dominate in mos t 



cases (Mart mez- Serrano et al."2009 Roskar et al.||201Q ). 
Sanchez-Blazquez et al. (2009) present cosmological sim- 
ulations of disk galaxies which show breaks in their ex- 
ponential profiles. While they find significant stellar mi- 
gration (60% of stars beyond the profile break formed 
in the inner disk), they attribute the profile break and 
corresponding increase in stellar age beyond the break to 
changes in the star formation rate. 

Secular stellar migration, caused by resonant inter- 
actions with disk structure ([Sellwood fc Binney 2002 



Roskar et al.| 2008a| [Schonrich fc Binney ||2009, ,Roskar| 
et al.||2011|), has attracted a lot of attention in recent 
years m part because it violates a key assumption that 
has been built into decades of chemical evolution models: 



that stars live out their lives where they are born (|Tins- 
Teyl[T975l [Matteucci fc Francois|ri989l ICarigi||19961 |Chi 



appini et al.||1997f Boissier fc Prantzos||1999 ). Migration 
has many consequences; for instance it naturally explains 
the flatness and large scatter in t he age-metallicity rela- 
tion at the solar neie; hborhood ( [Sellwood fc Binney|2002 



2 



Yoachim et al. 



Roskar et al. 2008a 



Schonrich fc Binney|[2009p . Stellar 

migration has even been proposed t o account for (a frac- 
tion of) the Milky Way thi ck disk (Schonrich & Binney 



2009[ [Loebman et al.|[2QTT). 



UBservationally constraining how often stellar migra- 
tion actually occurs in nature is an important next step. 
In external galaxies, This is ideally done using resolved 
stellar populations which are possible with the Hubble 
Space Telescope. The single-exponential g alaxy NGC 300 
shows no evidence of stellar migra tion ([Gogarten et al. 
" 2010 J. Howe ver, both NGC 4244 (Id e jQ5gir"al.| |!^007( 
Hoskar et aL| |2008b) and NGC 7793 (Had burn-Smith et 
aT 2011, submitted) show clear evidence of migration 
in the form of increasing mean stellar age beyond the 
break. Unfortunately, the number of disk galaxies suit- 
able for such resolv ed stellar po pulation studies is not 
very large. Instead, Bakos et al.| (12008) stack SPSS im- 
ages of galaxie s with similar profiles from the Pohlen 



Trujillo (2006) sample to measure radial color gradients. 
They find that galaxies with down-bending breaks have a 
g — r color minimum at the break radius. They take this 
as evidence that the stellar populations change across the 
break, with the region beyond the break being redder and 
older. 

In this paper, we aim to spectroscopically constrain 
the ages and metallicities of stars in the outer disks of 
many nearby galaxies to directly test if the stars are 
formed recently in situ or reach the outer disk through 
radial migration. To this end, we use the George and 
Cynthia Mitchell Spectrograph (formerly VIRUS-P)( |Hill 
|et al. 2008). With a large field of view, and large fibers, 
the Mitchell Spectrograph has proven to be an ideal 
instrument for observing nearby spiral galaxies (Blanc 
et al.|[2QQ9 ) and excellent at sp ectroscopically obser ving 



low surfac e brightness regions (iMurphy et al. 20111) . In 
Yoachim et al. (2010) we presented the analysis of one 
galaxy from our sample, NGC 6155, which showed an 
age upturn beyond the break. This paper presents our 
full sample. 

2. DATA 
2.1. Observations 

We obtained spatially resolved spectroscopy of our 
target galaxies using the Mitchell Spectrograph on the 
2.7m Harlan J. Smith telescope at McDonald Observa- 
tory. Observations were made during dark time in April, 
June, and August 2008 and March 2009. Identical instru- 
ment setups were used for all observations. The Mitchell 
Spectrograph with the VP-2 IFU bundle used here has 
a square array of 246 optical fibers which sample the 
1.9^x1.9' field-of-view with a 1/3 filling factor, with each 
fiber having a diameter of 4.23'^ The spectrograph im- 
ages the spectrum of each fiber on a 2048x2048 Fairchild 
Imaging CCD. On readout, the chip was binned by a fac- 
tor of two in the horizontal (spectral) dimension. 

For each galaxy, observations were taken at three 
unique dither positions to provide nearly complete spa- 
tial coverage. For the more extended galaxies, we also 
interleaved observations of blank sky. Individual galaxy 
exposures were 20 or 30 minutes with sky exposures of 
4-5 minutes. Twilight flat field frames, bias frames, and 
wavelength calibration Hg and Cd lamps were taken dur- 
ing the day. Flux standard stars were observed in twi- 
light. 



Our targets are listed in Table [T] We selected galaxies 
to match the instrument field of view and span a range of 
properties, including some with downward breaking pho- 
tometric profiles and so me with no break . Targets were 
primarily drawn from Pohlen & Trujillo (2006), wi th a 
few additional targets selected from HyperLeda ( |Paturel| 
eFar][2QQ3| ). 

2.2. Data Reduction 

We use the Mitchell Spectrograph software pipeline 
Vaccine (Adams et al. 2011) to perform basic reduction 
tasks. Vaccine subtracts the over-scan region and gener- 
ates and subtracts a bias frame. Bad pixels are masked. 
Twilight flats and arc lamps from each night are com- 
bined. The combined flat frame is used to trace the spa- 
tial profile of each fiber and correct for pixel-to-pixel vari- 
ation. The combined arc lamp is used to fit a wavelength 
solution for each fiber. We use a 5th order polynomial 
and find the wavelength solutions have a typical RMS of 
0.06-0.18 A 0.1 pixel). Vaccine then masks potential 
cosmic rays and collapses each fiber to a one dimensional 
spectrum. The fibers are adequately spaced, and our tar- 
gets are predominately low surface brightness, so we find 
making a fiber cross-talk correction unnecessary. 

Once the Vaccine routines are completed we make the 
following additional reductions using standard IRAF rou- 
tines and custom IDL scripts. The reduced twilight 
frames are used to generate a fiber-to-fiber flux correc- 
tion (similar to a longslit illumination c orrection). We 
use the b-spline procedure described in |Kelson| ( 2003 ) 
to generate background frames. For galaxy images, the 
bracketing sky observations are used to generate the sky 
frame while flux standard stars and smaller galaxies have 
enough empty fibers to generate a sky spectrum. Because 
we bin fibers together, imperfect sky-subtraction repre- 
sents a serious source of potential systematic errors. By 
taking sky observations between science exposures, we 
are able to construct background images that are nearly 
perfectly matched to the science frames, mimicking the 
success that has been achieved with nod- and- shuffle tech- 
niques. We have checked for over /under background sub- 
traction by cross-correlating binned spectra with the sky 
spectrum. If the sky was over (under) subtracted, the 
cross correlation would show a large maximum (mini- 
mum) at zero shift. We found no significant extrema in 
our cross-correlations. 

The spectra are then rectified to a common wavelength 
scale, and the flux photo-spectroscopic standards are 
used to flux calibrate the spectra and correct for atmo- 
spheric extinction. The spectra are then velocity shifted 
to the Local Standard of Rest. Images taken at iden- 
tical dither positions are averaged together with outlier 
rejection to eliminate any remaining cosmic rays. Many 
observations were taken in non-photometric conditions, 
so before combining frames we scale them to a common 
average flux level. 

Our final reduced data for each galaxy includes 738 
spectra spanning a 1.9^x1.9' field-of-view. The spectra 
have a wavelength range of ^3550-5580 with 2.2 A pix- 
els (with some fibers only reaching 5500 A due to the 
warping of the trace). The FWHM of reduced arc lamp 
images is 5.3 A, giving us a resolution of R ^ 700 — 1000. 
The fiber-to-fiber changes in resolution are less than 10%. 

The Mitchell Spectrograph utilizes a gimbal mount. 
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TABLE 1 
Galaxy Sample 



Name 


RA 


Dec 


Inclination 


Exposure Time 






Profile Type 




(2000) 


(2000) 


(degrees) 


degrees E of N 


(minutes) 


(km s-l) 


(km s-l) 


IC 1132 


15:40:07 


+20:40:50.1 


8 


140 


4555 


105.6 


II 


NGC 1058 


02:43:30 


+37:20:27.0 


18 


180 


746 


33.7 


I 


NGC 2684 


08:54:54 


+49:09:37.5 


22 


240 


2943 


101.4 


II 


NGC 2942 


09:39:08 


+34:00:22.6 


47 


180 


4412 


177.3 


II 


NGG 3888 


11:47:34 


+55:58:02.3 


47 


180 


2496 


203.1 


I 


NGG 4904 


13:00:59 


-00:01:38.8 


54 


180 


1025 


105.2 


II 


NGG 5624 


14:26:35 


+51:35:08.1 


56 


180 


2042 


66.9 


III 


NGG 6060 


16:05:52 


+21:29:06.0 


61 


280 


4503 


252.4 


II 


NGG 6155 


16:26:08 


+48:22:00.4 


34 


240 


2582 


109.8 


II 


NGG 6691 


18:39:12 


+55:38:29.5 


38 


180 


6130 


77.6 


II 


NGG 7437 


22:58:10 


+14:18:31.1 


23 


180 


2360 


151.9 


II 


UGG 04713 


09:00:20 


+52:29:39.3 


40 


180 


9127 


347.6 


II 



^ V and Vrot values taken from Hyperleda 

such that the instrument swings freely while the tele- 
scope slews. This mount keeps the gravity vector on the 
instrument constant through the night, resulting in excel- 
lent stability. In Figure [l] we show Lick index equivalent 
width measurements made from twilight frames. The 
fiber-to-fiber standard deviation of the Lick absorption 
equivalent widths is - 0.04 - 0.06 A, or 0.7-1.8%. This 
is excellent uniformity and stability, especially compared 
to earlier IFU studies where variations were of order 18% 
(e.g.. Figure 3 in |Ganda et al.||2007| . 

With our final redu ced data cubes, we used Pena lized 
Pixel-Fitting (pPXF) (|Cappellari fc Emsellem|2004|) and 



a|2004ft_ 



Gas AND A bsorption Line l^ltting (GANDALI^^) fSarzi' 
et al. [2006]) on each fiber to fit stellar velocities, stel- 
lar dispersions, and emission line velocities and fiuxes. 
For stellar templates we u se a suite of SSP spectra from 
Bruzual fc Chariot] ( [2QQ3| . We use GANDALF to mea- 
sure the [O IIJ doublet, [O III] doublet, and the higher 
order Balmer lines H/3, H7, H(5, He while masking the 
night skyline at 5577 A. The fit also includes a 4th de- 
gree multiplicative polynomial to correct for any dust 
or fiux calibration mis-match between the observations 
and templates. Because our data do not extend out to 
Ha, we cannot use the Balmer decrement to place tight 
constraints on the dust extinction. While pPXF fits the 
stellar velocity dispersion, our instrumental resolution is 
low enough (~ 140 km s~^) that we are unable to make 
meaningful velocity dispersion measurements beyond the 
central bulges of the most massive galaxies in the sample. 

The photometric fiux and gaseous velocity fields were 
used to fit a tilted rotating disk model to each galaxy. 
Once we have fixed the galaxy inclination and position 
angle, we median-filter the fiber surface brightness mea- 
surements and fit a broken-exponential profile (excluding 
the central 1-3 fibers which often have excess fiux from 
a bulge component). Our best-fitting profile parameters 
are listed in Table |2] The location of the profile break 
is well constrained, with uncertainties of l-2'^ The one 
exception is NGC 4904 where the inclination is poorly 
constrained, making the break radius fit unstable. The 
velocity fits were then used to shift each fiber to zero 
velocity. We assign a projected radius to each fiber, 
based on the best-fit galaxy inclination and position an- 
gle. Fibers with bright foreground/background objects 
were masked, and fibers were binned in elliptical annuli 
of 4 arcseconds. We then run the binned spectra through 



our star formation fitting procedure (described below). 

3. STAR FORMATION HISTORIES 

Stellar ages and metallicities have long been estimated 
using prominent spectral absorption features. The Lick 
indices have been u sed to go od effect on old elliptical sys- 
tems (e.g., Worthey ~t al.|1 994; Trag er et al.|2Q0Q] [Proc- 
'tor fc Sansom 2002 ; Sansom fc Northeast 2008). Unfortu- 
nately, measurements of Lick indices become difficult to 
interpret when the stellar population ca nnot be approx- 
imated as a s ingle-burst (S erra fc Trager|[2007 Yoachim 



fc Dalcanton 2008). When using individual equivalent 
widths (or broad band colors), one quickly runs into the 
problem of having more variables to fit than actual data 
points. 

Multiple authors have now attempted to model inte- 
grated galaxy sp ectra as the linear combination of SSP 
populations (e.g., Cid Fernandes et al.'2005^'Ocvirk et al 



2006; Mac Arthur et al.'2009; To jeiro et al. 2007; Koleva 
et al. 2009; Chihngarian et al. 2007| ISanchez et al.|201ip . 



See C hen et al., ( ,2010, ) for a comparison of many of these 
codes. 

A fundamental problem in fitting multiple Simple Stel- 
lar Populations (SSP) spectra to an observed galaxy is 
answering the question of how many spectra one needs? 
At low signal-to-noise in particular, the fraction of star 
formation in a particular age bin can become very un- 
certain. This can result in wildly different values for the 
average age of a population depending on the spectral 
range, signal-to- noise (SNR), and number of SSP models 
used in the decomposition. 

Inspired by the star formation histories derived from 
HST observations of resolved galaxies, we adopt model 
spectra that are continuous, rather than a combination of 
discrete SSPs. One important result of HST observations 
is that all g alaxies host old stellar populations at all radii 
(Wilhams et al ||2009alb| [Gogarten et al.||2010[ [Barker 
|et al.|i2007] ). This implies that any 'Sb ti decomposition 
that aoes not include old stars is probably ffawed. 

We continue the technique introduced in [Yoachim et al.| 
(2010), and explore the parameter space of exponentially 
rising and falling star formation histories. This technique 
should be expected to perform poorly in cases where 
the star formation histories are particularly bursty (e.g., 
dwarf galaxies or r egions of barred galaxies, Sanc hez^ 
[Blazquez et al.[ ( [2Qri, )), but we still expect to out-perform 
simple SSP decomposition. While hierarchical merging 
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Fig. 1. — Lick index equivalent widths measured from a twilight frame reduced identically to our science frames. The scatter in the 
values is very low and shows no overall spatial structure demonstrating the extreme stability of the instrument. Each point represents a 
single fiber. 



might lead to bursty star formation histories, we are 
primarily interested in the outskirts of spiral galaxies 
where observations of resolved stars have shown rela- 
tively smooth SFHs ( [Gogarten et al.||2QlQ ). 



3.1. Spectral Templates 



We use the stellar syntheses code Bruzual (2007) 
which updates Bruzual & Chariot ( 2003p to include ei- 
fects of TP-AGB stars. We have lliid the GALEXV 
package to generate a large suite of template spectra. 
The parameter space of the models is plotted in Fig- 
ure |2l The GALEXV code allows arbitrary SFHs, but 
includes only seven metallicity models. We have gener- 
ated 6 additional metallicities by interpolating between 
the GALEXV models. 

3.2. How Well Can Ages and Metallicities be 
Recovered? 

While there are numerous spectral synthesis codes, 
few have actually compared SSP decomposition to fit- 
ting more complicated SFHs. We have therefore r un a 
series of tests on the GANDALF ( |Sarzi et al.||2006[ ) rou- 
tine to find the optimal way to measure star formation 
histories. The simplest method is to find the best fit- 
ting simple stellar population (SSP). More complicated 
is to use a linear combination of multiple bursts of star 
formation (mSSP). Finally, one can use continuous star 
formation, usually parametrized by an exponential decay 
(rSP). 

We construct templat es using the GALEXV code 
( [Bruzual fc Charlot|2003 ). In particular, we constructed 



models with exponentially increasing or decreasing star 
formation, and models with constant star formation. We 
then ask what is the best way to measure a fiux- weighted 
average age — find an optimal combination of SSP mod- 
els or find the single best complex star formation history? 
We plot our templates in Figure [3] 

We construct artificial input spectra with known age 
and metallicity, re-binned to match our observations and 
with Gaussian noise added so that the spectra have a 
SNR per pixel of 50. For these tests, we use the wave- 
length range of 3500-5800 A and wavelength binning and 
resolution of 146 km s~^(2.2 A) and 5.3 A (the same as 
our data cubes). 

In Figure [4j we show our results of using GANDALF 
to fit an SSP spectrum with a linear combination of mul- 
tiple SSP (mSSP) template spectra. We point out that 
this is an idealized fitting scenario. The input spectra 
and template spectra are from the same source, we have 
not added any emission lines and there is no template 
mismatch. For old ages and high metallicities, the fits 
are excellent. However, the fits diverge when the metal- 
licity is sub-solar or the fiux- weighted age is less than 100 
Myr. This is a fairly damning experiment for trying to 
decompose SFHs with mSSP as it implies just the pres- 
ence of Gaussian noise is enough to trick the code into 
using more templates than are actually necessary. In- 
specting the P^i" degree-of-freedom for these fits, we 
find all the fits fell between 0.9 and 1.1. There is no par- 
ticular correlation between the values and how well 
the fits matched the input age and metallicity. 

Next, we test how well mSSPs can fit rSP spectra. Our 
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Fig. 2. — Parameter space covered by our spectral synthesis models. We span a large range of metallicities and include very old stellar 
populations as well as populations dominated by recent formation. 
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Fig. 3. — Our template spectra, (a) Models with constant star formation rates, (b) Model spectra from declining star formation rates, 
(c) Model spectra generated from increasing star formation rates. 



results are plotted in Figure [5| where we fit increasing, 
decreasing, and constant star formation histories. As 
before, the mSSP work well for old metal-rich popula- 
tions. The fits are biased to returning young ages and 
high metallicities for the majority of input spectra. This 
is caused by GANDALF's propensity to often place a 
small amount of weight on the oldest stars. Again, the 
per degree-of- freedom fell between 0.88 and 1.13 for 
these fits. 

The results of decomposing an observed spectrum into 
multiple SSP models is also very wavelength dependent. 
If we limit the SSP decomposition to the region redward 
of 4040 A, we find the resulting average age can be two or 
more Gyr different from fitting the full spectrum. This 
is a well known effect where the Ca II H and K lines ar e 
v ery sensitive to young populatio ns ( Smith et al.||2009 ). 



Marmol-Queralto et al. (2011) use a similar analysis 
and also investigate varying the SNR. They also find 
that decomposing with SSP models is difficult, with flux- 
weighted age errors of ~ 50% for old stars and metallicity 
errors of ~ 100%. 

We must emphasize that this is not a failure on the 
part of the GANDALF code-if anything we are guilty of 
not using GANDALF the way it was intended. GAN- 
DALF was written primarily to make better emission 
line measurements while correcting for underlying stellar 
absorption features. 

Figure [6] shows the results of fitting extended star for- 
mation histories with the best single template. Unlike 
fitting a combination of SSPs, this technique is very ro- 
bust against varying SNR and wavelength ranges. The 
major disadvantage is that the fits should be expected to 
return poor fits for galaxies which have bursty star for- 



mation histories. We also do not include any metallicity 
evolution, assuming all the stars have a similar metallic- 
ity. While this sounds like a rather drastic assumption. 



Williams et al. ([2009b) use resolved stars in the outskirts 
of M81 and tina that over the entire star formation his- 
tory of their observed region the average metallicity only 
varies by ~ 0.5 dex. Figure [7] shows results from a simu- 
lation where we added simple chemical evolution to the 
spectra to be fit. As expected, our fitting procedure does 
a poor job matching the metallicity, however, we still do 
an excellent job recovering the ages. Figure [7| also shows 
the results of fitting star formation histories that are not 
perfect exponentials. In these cases, we added varying 
amounts of noise to smooth exponential star formation 
histories with bursts of one- hundredth, one-tenth, and 
one-third the smooth SFR. Our fitting method returned 
excellent results, failing only when there was a recent 
strong burst of star formation. 

Based on our experience fitting HST derived star for- 
mation histories, we find our technique has systematic 
error of ~0.6 Gyr, a nd that the Ax^ = 20 returns a rea- 
sonable uncertainty ( Yoachim et al.||2010 ). Examples of 
our preferred fitting procedure are shown in Figure [8) 
Because our residuals are dominated by systematic dif- 
ferences (rather than Gaussian noise), we cannot legiti- 
mately use formal uncertainties on the fitted param- 
eters. 

Figure |9] shows examples of our binned spectra for two 
galaxies. In particular, we show how GANDALF can be 
used to remove emission lines to recover the underlying 
stellar features. Based on our previous experience, we 
only attempt to fit spectra where the SNR per pixel is 
greater than 40. 
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SSP with mSSP SSP with mSSP 




Input age (yrs, flux weighted) Input [Z/H] (flux weighted) 

Fig. 4. — Results of fitting single SSP spectra with mSSP models. 





Fig. 6. — Fitting extended star formation histories with the best single extended SFH template. Solid lines show the perfect-fit case. 
Unlike fitting with a linear combination of SSP spectra, we do an excellent job of recovering the input parameters. 
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Fig. 7. — Left: Results from fitting star formation histories which have experienced significant chemical evolution. While we fail to 
recover the correct metallicities, our fitted ages are still quite robust. Right: Fitting star formation histories which are not perfectly smooth 
exponentials. Our fits only fail when the recent star formation is very bursty. 
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Fig. 8. — Example of our stellar population fitting results. Top panels show the surface generated using GANDALF and a set of 
exponential star formation histories. Contours mark the Ax^ = 2.3,20 (thick), and 50 levels. Model points are marked with black circles 
and the best fit models are marked with stars. The lower panels show the best fit model (stellar spectrum plus Gaussian emission lines) in 
red along with the raw spectrum and residuals in black. 



8 



Yoachim et al. 




Fig. 9. — Examples of our spectra after binning fibers. Top panels show the binned spectra while bottom panels show the results after 
subtracting off the best-fit emission lines. 
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4. ANALYSIS 

We present the galaxies ordered by their maximum ro- 
tational velocity. 

4.1. Notes on Individual Galaxies 

4.1.1. NGC 1058 

The galaxy nearly fills the field of view and is close to 
face-on, resulting in limited rotation curve information. 
We do not detect any break in the surface brightness 
profile. The stellar age profile shows a smooth transition 
from an old center to a young star forming outer region. 
The age profile levels off and is constant in the outer 
region of the galaxy. There is a surprising amount of 
spiral structure for such a low mass system. 

4.1.2. NGC 5624 

The only Type III (up-bending) profile in our sam- 
ple. Our stellar population fits show this system has 
recently had a central burst of star formation, while the 
outer region is significantly older. The H/3 emission is 
very concentrated in the inner region of the galaxy. Our 
observations were not deep enough to probe the region 
beyond the profile break. Similar central star bursts have 
been found in other dwarf g alaxies with up-bending pro- 
files ( Herrmann et al.||2Qll ). Anti-truncation via a cen- 
tral burst is al s o in a greement with the simulations of 
[Younger et aLj (2007) who use minor- mergers to drive 
gas to the center of galaxies to create an anti-truncated 
surface brightness profile. 

4.1.3. NGC 6691 

The best-fit age does change slightly beyond the break 
radius, however, we do not consider this strong evidence 
for migration as the D4000 indicator shows little change 
across the break. The stellar populations beyond the 
break are also still best fit by an increasing star forma- 
tion rate, implying that the profile break does not cor- 
respond to a star formation threshold. This could be a 
case where a recent burst of star formation in the outer 
disk is masking the dominant older stellar population. 
Such bursty star formation histories will be poorly fit by 
our technique. 



4.1.4. NGC 



[Pohlen fc Trujillol ([2006) fit the break for this galaxy 
at 24'' and classify it as a Type H-CT (classical trun- 
cation). We find a similar break radius, although it is 
not particularly pronounced. There is an increase in the 
stellar age around the break radius. Interestingly, we 
find the age minimum is actually well inside the break 
radius at around lb" . The few stellar velocity points we 
have beyond this radius still show rotation. There is also 
a strong metallicity gradient, with the outskirts being 
metal-poor compared to the interior. 

4.1.5. NGC 4904 

Our ellipse fitting for the photometry of this galaxy 
is not stable, with the center and position angle of the 
galaxy being difficult to pin down. The kinematic po- 
sition angle is also noticeably different from the photo- 
metric major axis. Spectroscopically, it looks like active 
star formation across the disk, and also outside the break 



radius. Pohlen & Trujillo (2006) find the break radius at 
39^' in r and 104" in g and suggest the break is due to 
the outer Lindblad resona nce. Our best-fit break ra dius 
is around 49'', between the [Pohlen fc Trujillo[ ([2006[) val- 
ues. The galaxy has a strong bar, and visually looks like 
the break is near the expected OLR radius. 

4.1.6. IC 1132 

This galaxy shows a very strong profile break, but no 
significant change in the stellar population across the 
break. Many fibers show H/3 emission beyond the break 
radius as well, suggesting robust star formation beyond 
the break radius. 

4.1.7. NGC 6155 

A strong inner star forming region creates a minimum 
in the stellar age. There is a clear increase in the ^ie\- 
lar ages beyond the profi le break as originally reported 
in Yoachim et al.[ (2010). We also detect a decrease in 
the stellar metallicity beyond the break radius, consistent 
with the migration hypothesis. There is little Hy ^ emis- 
sion or UV-emiss ion beyond the break radius. [Pohlen 



Trujillo 



20061 ) fit the break at 34'', while our fit is 
a comparable 36 . Again, the stellar kinematics remain 
disk-like beyond the break radius. 

4.1.8. NGC 1431 

The stellar ages increase beyond the profile break, but 
only one point is fit with a decreasing star formation 
history. Similar to NGC 2684, the actual age minimum 
actually occurs inside the break radius. Emission from 
H/3 does get sparse beyond the break radius. This galaxy 
is liste d as a classical truncation in Pohlen fc Trujillo] 
(2006) where they find a similar break radius to our fit. 



4.1.9. NGC 2942 

The break is not particularly well defined as bright 
spiral arms extend well beyond our best-fit break radius. 
Spectroscopy failed to reach adequate signal-to-noise be- 
yond the break. The SDSS and GALEX imaging shows 
extensive star formation in the spiral arms extending be- 
yond the break. We also detect H/3 emission beyond the 
break. 

4.1.10. NGC 3888 

The surface brightness profile of NGC 3888 is well de- 
scribes by a singl e exponential, as also found in Pohlen 
fc Trujillo (2006). The stellar age profile shows a pro- 
nounced older population around 30". 

4.1.11. NGC 6060 

We failed to reach adequate SNR to measure a spectro- 
scopic age beyond the break. The GALEX and H/3 im- 
ages do suggest that star formation truncates at the sur- 
face brightness profile break. This is a high mass galaxy 
(vrot=252 km s~^) with a very smoothly declining stellar 
age profile. The best fit SFHs are all decreasing, implying 
the gas supply for this galaxy is be ing exhausted. This 
galaxy was not in 'Pohlen fc Trujillo ( 2006 ) . While we 
detect a break radius, we did not have adequate signal 
to measure stellar populations beyond the break radius. 
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Fig. 10. — NGC 1058. Upper left: SDSS or DSS image of the target. Upper Middle: V-band surface brightness measured from Mitcheh 
Spectrograph spectra. Fibers below a signal-to-noise threshold have been masked along with those containing bright foreground stars. 
Upper Right: Near UV image from GALEX (if available). Middle Left: Stehar velocity field measured with pPXF. Middle Middle: H/3 
emission line fiux after correcting for stellar absorption. Middle Right: Radial surface brightness profile measured from fiber magnitudes. 
A yellow line marks a running median, the solid green line shows the best fit broken exponential. The broken green lines shows single 
exponentials fit inside and outside the break radius and extrapolated to the full radial range. Lower Left: best-fit mass-weighted stellar 
ages from binned spectra. Red points are used for models that converged on exponentially decreasing SFHs, while blue points are used for 
increasing SFHs. We also show the age sensitive D4000 index for comparison (The D4000 index has not been corrected for dust effects). 
Lower Middle: The best fit values of the star formation timescale 1/r, along with the emission-corrected Balmer line absorption equivalent 
widths. Lower Right: Best-fit fiux-weighted average metallicities. In most figures, a red line marks the surface brightness break radius. 
For galaxies without a break radius, we mark the iiy = 23 mag/D^^ contour with a black line. Error bars are derived from the Ax^ = 20 
contour (with a systematic error fioor of 0.6 Gyr for the average age plots). 



4.1.12. UGC 04713 

There is a slight upturn in stehar ages before the profile 
break, suggestive of possible stellar migration. The H/3 
emission and UV flux are truncated well before the profile 
break radius. Like NGC 6060, this is a high mass system 
(vrot=347 km s~^) where all the SFHs are decreasing, 
implying very low levels of star formation. Again, we 
did not reach adequate signal-to-noise to measure stellar 
populations beyond the break. 

5. DISCUSSION 

Figures P^OpTI show our observations and best-fit stel- 
lar ages and metallicities. We have ordered the galaxies 
by rotation velocity. Three of the galaxies, listed in Ta- 
ble |3j show indications of stellar migration with the stars 
beyond the profile break being 1-2 Gyr older than those 
inside the break. The star formation history for these 
outer disks are also best fit by declining rates, while the 
interiors show signs of active star formation and are best 
fit with increasing rates of star formation. In all of these 
galaxies, both the D4000 index and the Balmer absorp- 



tion lines show significant changes indicating the out- 
skirts of the galaxies are older than the region inside the 
break radius. 

Exploring low surface brightness levels, there is always 
the concern that the observations could become domi- 
nated by a stellar halo population. For the few galaxies 
where we have adequate SNR to measure stellar kine- 
matics beyond the break radius (NGC 2684 and NGC 
6155), the velocity field remains disk-like, suggesting we 
have not become dominated by a halo or stream stellar 
population. 

While we find several examples where the surface 
brightness profile break corresponds with a change in 
stellar population, the actual minimum is rarely at the 
break radius. Rather, the stellar age minimum is often 
located well inside the break radius. In most cases, it 
appears the region inside the break radius is hosting the 
most active star formation, thus driving down the aver- 
age age compared to the SFR near the profile break. It 
should be emphasized that these results are in systems 
with Vrot ^ 100 km s~^, and may not apply to higher 
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Fig. 12. NGC 6691. Same as Figure [Tol 
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Fig. 14. NGC 4904. Same as Figure [To] 



Resolved Star Formation Histories 




Fig. 16. NGC 6155. Same as Figure [To] 
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Fig. 18. NGC 2942. Same as Figure [To] 
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Fig. 19. — NGC 3888. Same as Figure [Tol The black ellipse marks the muy = 23 contour. 
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Fig. 20. NGC 6060. Same as Figure [To] 
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Fig. 21. UGC 04713. Same as Figure [To] 



TABLE 2 
Fitted Structural Parameters 



Name 




II 


II 


Pouter 
II 


IC 1132 


21.1 


17.0 


26.6 


7.4 


NGC 1058 


20.5 


25.8 






NGC 2684 


20.1 


9.0 


22.3 


6.7 


NGC 2684^ 


19.6 


10.6 


24 


5.6 


NGC 2942 


20.8 


17.3 


45.5 


10.8 


NGC 3888 


19.0 


11.0 






NGC 3888^ 


18.7 


10.1 






NGC 4904 


20.3 


23.2 


48.8 


11.1 


NGC 4904^ 


20.2 


27.6 


39 


12.5 


NGC 5624 


19.8 


9.2 


39.3 


21.0 


NGC 5624^ 


20.1 


9.9 


48 


19.8 


NGC 6060 


20.1 


21.3 


56.6 


12.3 


NGC 6155 


19.7 


11.6 


35.6 


7.8 


NGC 6I55I 


19.4 


12.2 


34 


8.1 


NGC 6691 


20.6 


16.1 


36.5 


8.6 


NGC 7437 


21.6 


24.2 


43.6 


9.7 


NGC 7437I 


21.3 


26.4 


43 


13.9 


UGC 04713 


20.2 


13.3 


51.1 


9.5 



r-band fits from |Pohlen Trujillo] ( |2006| 
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TABLE 3 

Properties of Migration Signature Galaxies 



dominate. 



6. CONCLUSIONS 



Name 


Break Strength 


A age 


A Z 


Vrot 




hinner / houter 


(Gyr) 


(Dex) 


(km s-1) 


NGC 2684 


1.40 


1.1 


0.00 


101 


NGC 6155 


1.46 


2.4 


-0.46 


110 


NGC 7437 


2.46 


1.1 


-0.11 


152 




TABLE 4 






Properties of No Migration Signature Galaxies 


Name 


Break Strength 


A age 


A Z 


Vrot 




hinner / houter 


(Gyr) 


(Dex) 


(km s-1) 


IC1132 


2.31 


0.2 


0.0 


106 


NGC4904 


2.10 


0.5 


-0.18 


105 


NGC6691 


1.57 


-0.1 


-0.29 


78 



mass systems. Since the spectra can be dominated by 
bright current star formation, we could be biased to mea- 
suring younger ages by strong current star formation. Al- 
ternatively, star formation may not proceed in a simple 
inside-out fashion, and the complex star formation his- 
tories of real galaxies can move the stellar age-minimum 
away from the break radius. 

While three galaxies in our sample show signs of stellar 
migration, three of the galaxies (listed in Table [4|) show 
little or no sign of the outer disk being formed through 
stellar migration despite having some of the strongest 
profile breaks. While there is no strong upturn in the 
age profiles of these galaxies, we also note that the ages 
do not seem to be decreasing either. 

Unfortunately, we failed to reach adequate signal-to- 
noise in higher mass systems, so we cannot comment on 
how efficient stellar migration is as a function of galaxy 
mass. Interestingly, the galaxies with the weakest pro- 
file breaks show the strongest evidence for stellar mi- 
gration, while galaxies with stronger profile breaks have 
outer regions dominated by younger stars. This could be 
consistent with the migration hypothesis, as the galaxies 
that have experienced the most migration should have 
the weakest profile breaks. 

Our metallicity measurement is rather crude since we 
do not include a chemical enrichment model with the star 
formation histories. It is notable that both NGC 2684 
and NGC 6155 show both an age increase correlated with 
a metallicity decrease. However, NGC 4904 and NGC 
6691 also show a decrease in metall icity past the break 
radius but no change in stellar age. |Vlajic et al.| (|2011|) 
look at the outer disk of NGC 300 which is well described 
by a single exponential with resolved stars and find an 
upturn in the metallicity gradient. 

Overall, we find a very wide range of stellar populations 
beyond profile breaks, with some outer disks dominated 
by old stars while others host vigorous star formation. 
It would thus appear profile breaks in spiral galaxies will 
require multiple formation mechanisms. For example, ra- 
dial migration may be the primary driver of stellar popu- 
lation trends across the profile break for isolated galaxies, 
while recent infall of gas may drive star formation in the 
outskirts of a galaxy that generates a break. Cosmologi- 
cal simulations should help determine which mechanisms 



We spectroscopically confirm that stellar population 
ages increase near spiral galaxy surface brightness profile 
breaks. However, we also find several examples where 
star formation is clearly active on both sides of the break 
radius and the outer disk is dominated by young stars. 

While Bakos et al.^ (2008) found that, on average, the 
break radius corresponds to a change in stellar age, we see 
that in individual cases the stellar population can grow 
older or remain constant. Our results suggest that the 
behavior of profile breaks is bi-modal, with some galaxies 
hosting a change in stellar populations, while others have 
similar stellar popu lations throughout the disk. This is 
consistent with the 



Bakos et al. 



(2008) result, but sug- 



gests that the population change they detect has been 
diluted by galaxies where there is no change in the stel- 
lar populations beyond the breaks. 

In the future, we would like to extend our SFH fitting 
procedure to fit multiple epochs of star formation rather 
than being limited to continuous functions. We would 
also like to expand to further constrain the star forma- 
tion histories with broad-band colors, possibly combining 
optical spectra with IR and UV observations. 

There is obviously a wealth of data present in IFU 
observations. While we have focused on stellar popu- 
lations, our data also allow the construction of velocity 
maps and measurements of emission line strengths. We 
plan to release our re duced data-cubes along with the 
larger VENGA survey dBlanc et al.||2010| ). The data are 
currently available on tne author's web site^. 
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